Two and One-dimensional Honeycomb Structure of Boron Nitride 
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This paper presents a systematic study of two and one dimensional honeycomb structure of boron 
nitride (BN) using first-principles plane wave method. Two-dimensional (2D) graphene like BN 
is a wide band gap semiconductor with ionic bonding. Phonon dispersion curves demonstrate the 
stability of 2D BN flakes. Quasi ID armchair BN nanoribbon are nonmagnetic semiconductors 
with edge states. Upon passivation of B and N with hydrogen atoms these edge states disappear 
and band gap increases. Bare zigzag BN nanoribbons are metallic, but become a ferromagnetic 
semiconductor when their both edges are passivated with hydrogen. However, their magnetic ground 
state, electronic band structure and band gap are found to be strongly dependent on whether B- or N- 
edge of the ribbon is saturated with hydrogen. Vacancy defects in armchair and zigzag nanoribbons 
affects also magnetic state and electronic structure. In order to reveal dimensionality effects these 
properties are contrasted with those of various 3D BN crystals and ID BN atomic chain. 

PACS numbers: 73.22.-f, 75.75.4-a, 63.22.-m 
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I. INTRODUCTION 

Synthesis of a single atomic plane of graphite, i.e. 
Graphene with covalently bonded honeycomb lattice has 
been a breakthrough for several reasonsii^i^. Firstly, 
electrons behaving as if massless Dirac Fermions have 
made the observation of several relativistic effects pos- 
sible. Secondly, stable graphene has disproved previ- 
ous theories, which were concluded that two-dimensional 
structures cannot be stable. Graphene displaying excep- 
tional properties, such as high mobility even at room 
temperature, ambipolar effect, Klein tunelling, anoma- 
lous quantum hall effect etc. seems to offer novel appli- 
cations in various fields^. Not only 2D graphene, but 
also its quasi ID forms, such as armchair and zigzag 
nanoribbons have shown novel electronic and magnetic 
properties^ii, which can lead to important applications 
in nanotechnology. As a result, 2D honeycomb struc- 
tures derived from Group IV elements and Group III-V 
and II- VI compounds are currently generating significant 
interest owing to their unique properties. 

Boron-Nitride (BN) in ionic honeycomb lattice which 
is the Group III-V analogue of graphene have also 
been produced having desired insulator characteristics^. 
Nanosheet o^'^" , nanoconcs^^, nanotubes^^, nanohorns^'^, 
nanorodsi^ and nanowiresi^ of BN have already been 
synthesized and these systems might hold promise for 
novel technological applications. Among all these differ- 
ent structures, BN nanoribbons, where the charge car- 
riers are confined in two dimension and free to move in 
third direction, are particularly important due to their 
well defined geometry and possible ease of manipulation. 

BN nanoribbons posses different electronic and mag- 
netic properties depending on their size and edge ter- 
mination. Recently, the variation of band gaps of BN 
nanoribbons with their widths and Stark effect due to 
applied electric field have been studie d^^i^^ . Recently 
the magnetic properties of zigzag BN nanoribbons have 



been investigated^^. Half- metallic properties have been 
revealed from these studies which might be important for 
spintronic applications. Production of graphene nanorib- 
bons as small as 10 nm in width has been achieve d^^i^° 
and similar techniques are expected to be developed for 
BN nanoribbons. 



A thorough understanding of 2D BN honeycomb struc- 
ture and their various nanoribbons is important for fur- 
ther study of this graphene like compounds. BN by it- 
self provides with very interesting chemical and physical 
properties, which may lead to important applications. In 
this paper, we present a detailed ab-initio study of elec- 
tronic, magnetic and elastic properties of 2D (graphene 
like) BN and bare and hydrogen passivated, quasi ID 
BN nanoribbons (BNNRs). We also investigated the ef- 
fect of the vacancy defects on these properties. To reveal 
the dimensionality effects we include also a short discus- 
sion regarding 3D BN bulk crystals and ID BN atomic 
chains. We found that 2D BN is a nonmagnetic, wide 
band gap semiconductor. The ionic bonding due to sig- 
nificant amount of charge from B to N atom opens a 
gap and hence dominates electronic structure. Calcu- 
lated phonon dispersion curves provide a clear evidence 
that 2D BN flakes is stable. The armchair and zigzag 
nanoribbons of BN display even more interesting elec- 
tronic and magnetic properties. Bare and hydrogen pas- 
sivated armchair BN nanoribbons (A-BNNR) are non- 
magnetic wide band gap semiconductor. The value of 
band gap of A-BNNR having width w > 10 A is practi- 
cally independent from the width of nanoribbons. While 
the bare zigzag BN nanoribbons (Z-BNNR) are ferro- 
magnetic metal, they become nonmagnetic semiconduc- 
tor upon the passivation of both edges. We found that 
2D BN and its nanoribbons have properties, which are 
complementary to graphene. 
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II. MODEL AND METHODOLOGY 

We have performed first-principles plane wave calcu- 
lations within density functional theory (DFT) using 
PAW potentials-i. The exchange correlation potential 
has been approximated by generalized gradient approx- 
imation (GGA) using PW91^^ functional both for spin- 
polarized and spin-unpolarized cases. All structures have 
been treated within supercell geometry using the peri- 
odic boundary conditions. A plane-wave basis set with 
kinetic energy cutoff of 500 eV has been used. In the 
self-consistent potential and total energy calculations the 
Brillouin zone (BZ) is sampled by special k-points. The 
numbers of these k-points are (15x15x15) for bulk BN, 
(25x25x1) for 2D BN and (25x1x1) for nanoribbons, re- 
spectively, and are scaled according to the size of su- 
perlattices. All atomic positions and lattice constants 
are optimized by using the conjugate gradient method, 
where the total energy and atomic forces are minimized. 
The convergence for energy is chosen as 10^^ eV between 
two steps, and the maximum Hellmann-Feynman forces 
acting on each atom is less than 0.02 eV /A upon ionic 
relaxation. A large spacing (~ 10 A) between monolay- 
ers has been taken to prevent interactions between them. 
The pseudopotentials having 3 and 5 valence electrons 
for the B (B: 28^ 2p^) and N ions (N: 25^ 2p^) were 
used. Numerical calculations have been performed by 
using VASP package22i2i. The phonon dispersion curves 
are calculated within density functional perturbation the- 
ory (DFPT) using plane wave methods as implemented 
in PWSCF softwares^. 



III. 3D BN CRYSTALS AND ID ATOMIC 
CHAIN 

In this section, we present our theoretical calculations 
on 3D bulk BN crystals and truly ID BN atomic chain. 
Earlier these 3D bulk crystala^^i^i^^iSai^ and ID atomic 
chains^ have been studied theoretically by using differ- 
ent methods. Our purpose in including these crystals 
of BN in different dimensionalities is to contrast their 
properties with those of 2D and quasi ID honeycomb 
structures of BN and also reveal dimensionality effects. 

A. 3D Bulk BN crystals 

Three dimensional bulk crystals include hexagonal- 
layered BN (h-BN), wurtzite BN (wz-BN) and zincblende 



BN (zb-BN) structures. Their atomic configurations and 
primitive unit cells are described in Fig. [1] By using the 
expression, 

Ec = E[BN] - E[B] - E[N] (1) 

where E[BN] is the total energy per B-N pair of the op- 
timized structure of BN crystal; E[B] and E[N] are the 
total energies of free B and N atoms; we calculated the 
equilibrium cohesive energies of h-BN, wz-BN and zb-BN 
crystals as -17.65, -17.45 and -17.49 eV per B-N pair, re- 
spectively. Accordingly, h-BN, which is the analogue of 
graphite, is the most energetic bulk structure. On the 
other hand, the cubic BN structure is known to be the 
second hardest material of all. 

The lattice constants of the optimized structures in 
equilibrium are a = 2.511 A, c/a — 2.66 and the dis- 
tance between the nearest B and N atoms is d = 1.450 
Afor h-BN layered crystal. For wz-BN, optimized val- 
ues of a, c/a and d are calculated to be 2.542 A, 1.64, 
and 1.561 A, respectively. The zincblende structure has 
lattice constant a — 2.561 A and d = 1.568 A. All 
our results related with the structural parameters are 
in good agreement with the experimental and theoretical 
values^i^LSSiiSSi^ within the average error of ^ 1%. 

The calculated electronic band structure, total and 
partial (or orbital projected) density of states (DOS) of 
3D crystals are presented in Fig. [TJ These h-BN, wz-BN 
and zb-BN crystals are indirect band gap semiconductors 
with calculated band gaps being i?G=4.47, 5.72, and 4.50 
eV, respectively. The calculated values of Eg differ from 
the earlier ones depending on the method usedS. For 
h-BN having 2D BN atomic layers in the (x,y)-plane. 
The band structure is composed from the band struc- 
tures of these individual atomic layers with hexagonal 
symmetry, which are slightly split due to weak coupling 
between them. Highest valence band has N-p^ character; 
the states of the lowest conduction band is formed from 
B-pz orbitals (the z direction corresponds to "c" in Fig. [T] 
). Overall features of the total density of states (TDOS) 
are similar for three 3D crystal structures. Valence band 
consists of two parts separated by a wide intra band gap. 
The lower part at -20 eV is projected mainly to N-s 
and partly to N-p and B-s orbitals. The upper part is 
due to mainly N-p and partly B-p orbitals and has sim- 
ilarities in both zb-BN and wz-BN crystals. As for the 
lower part of the conduction band it is derived mainly 
from B-p orbitals. The differences of three 3D crystals 
are pronounced in the lower part of the conduction band. 



We calculate the amount of charge on constituent B 
and N atoms in 3D crystals by performing the Lowdin'^^ 
analysis in terms of the projection of plane-waves into 
atomic orbitals. By subtracting the valencies of free B 



and N atoms from the calculated charge values on the 
same atoms in 3D crystals we obtain the charge transfer, 
AQ. The calculated values of AQ for h-BN, wz-BN and 
zb-BN are 0.416, 0.342, 0.334 electrons, respectively The 
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FIG. 1: (Color online) Optimized atomic structure, energy bands, total (TDOS) and orbital projected density of states (PDOS) 
of various 3D crystals of BN. (a) Hexagonal (h-BN) whose B(N) atoms are on top of the N(B) atoms in the consecutive layer; (b) 
wurtzite (wz-BN); and (c) zincblende (zb-BN) crystals. Dark-green and light-gray balls represent B and N atoms, respectively. 
The band gaps between conduction and valence bands are highlighted. The orbital character of states are indicated for the 
conduction and valence band edges. The zero of energy is set to the Fermi energy Ef. All structures are fully optimized. 



fact that AQ of zb-BN and wz-BN have almost equal val- 
ues, but AQ of h-BN crystal is significantly larger related 
to the shorter B-N bond length in h-BN crystal. 



B. ID BN Atomic chain 

BN forms stable segments of linear atomic chain^ 
like carbon^^. This situation is in contrast to second 
and third row elements (such as Si and Ge) and III- 
V compounds and metals (such as Al, Au etc) which 
can form stable zigzag chain structures instead of lin- 
ear chain structures. Our results on optimized chain 



structure yield the cohesive energy £'c'=16.04 eV per 
B-N pair, the B-N distance d =1.307 A, the indirect 
band gap Eg =3.99 eV and charge transfer from B to N, 
AQ =0.511 electrons. Hence the double bond between B 
and N is ionic. 



IV. 2D HONEYCOMB STRUCTURE OF BN 

Having discussed the overall structural and elastic 
properties of 3D and ID BN, we now consider 2D BN 
with hexagonal symmetry. The atomic structure of 2D 
BN is similar to the honeycomb structure of graphene. 
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except that the constituent atoms of the former are from 
III and V columns of the Periodic Table. Normally, the 
bond between nearest B and N atoms is formed from the 
bonding combination of B-sp^ and N-sp^ orbitals. How- 
ever, owing to the electronegativity difference between B 
and N atoms electrons are transferred from B to N. As 
a result, in contrast to purely covalent bond in graphene 
the bonding between B and N gains an ionic character. 
The charge transfer from B to N dominates several prop- 
erties of 2D BN including the opening of the band gap. 
In this respect the BN honeycomb structure is comple- 
mentary to graphene. 



A. Charge density analysis and electronic structure 

The atomic structure, atomic charge, charge transfer 
from B to N and the electronic structure of 2D BN are 
presented in Fig. O Contour plots of total charge indi- 
cates high density around N atoms. The difference charge 
density is calculated by subtracting charge densities of 
free B and N atoms from the charge density of 2D BN, 
i.e. A/9 = pbn — Pb — Pn- High density contour plots 
around N atoms protruding towards the B-N bonds in- 
dicate charge transfer from B to N atoms. This way the 
B-N bonds achieve an ionic character. The amount of 
transfer of charge is calculated by Lowdin analysis to be 
A(3=0.429 electrons. Interestingly, AQ is slightly larger 
than that calculated for h-BN, but significant larger than 
those calculated for wz-BN and zb-BN crystals. 

2D BN is a semiconductor. Calculated electronic en- 
ergy bands are similar to those calculated for h-BN crys- 
tal. The TT- and tt*- bands of graphene which cross at 
the K- and K*-points of the BZ open a gap in 2D BN 
as a bonding and antibonding combination of N-p^ and 
B-P2 orbitals. The contribution of N-p^ is pronounced 
for the filled band at the edge of valence band. The cal- 
culated band gap is indirect and Eq =4.64 eV. TDOS 
and partial density of states show also similarity to those 
of h-BN layered crystal presented in Fig. [TJ 



B. Phonon spectrum 

Even if the structure optimization resulting in the hon- 
eycomb structure in Fig. [2] can be taken as an indication 
for the stability of 2D BN, calculation of phonon dis- 
persion curves through the diagonalization of dynamical 
matrix provides a more stringent test for stability. One 
of acoustical branches for F to X curves taking negative 
value even at a small region of BZ indicates the insta- 
bility of the structure. There have been a number of 
experimental® and theoretical studies of phonon spec- 
trum of 2D21 and 3D honeycomb B1<!M^S^^^'^^. Here, 
the phonon dispersion curves of h-BN, 2D BN and ID 
BN chain and density of states together with the infrared 
(IR) and Raman (R) active modes of 2D BN and h-BN 
at F-point have been calculated by using density func- 
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FIG. 2: (Color online) (a) Primitive unit cell of the hon- 
eycomb structure of 2D BN together with Bravais lattice 
vectors. Calculated total charge density pBN and difference 
charge density Ap, are also shown in the same panel, (b) 
Calculated electronic structure of 2D BN honeycomb crys- 
tal together with total, TDOS and partial density of states, 
PDOS on B and N atoms. The orbital character of the states 
are also indicated. 



tional perturbation theory (DFPT) as implemented in 
PWSCF software^. For the DFPT phonon calculation 
of bulk h-BN, we used a four atom primitive cell, which 
yield 12 phonon branches at the center of BZ in Fig. 
[3] (a). The symmetry point group is calculated as Dg/i 
(space group P6/mmm). The irreducible representations 
at F is 2 E2g-f 2 B2g-h2 A2„-h2 Ei„. While the modes 
Ei„ and E2g are doubly degenerate, B2g and A2ti are 
non degenerate. The modes Ei„ and A2ti are IR active, 
the E2g is Raman active. B20 is an inactive mode. Our 
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results are in agreement with previously calculated and 
experimental data, but differ slightly from those of Ser- 
rano et al.iS. While present GGA calculations predict 
B2g mode as an inactive mode, LDA calculations by Ser- 
rano et al. found Big as an inactive mode. Most of the 
phonon bands of h-BN are degenerate. This indicates 
that the coupling between BN layers in h-BN is weak. 
However, it is well known that the BN is polar mate- 
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rial with long range dipole-dipole interaction. This gives 
rise to the splitting between longitudinal optical (LO) 
and transverse optical (TO) mode at F point. The low- 
est transverse acoustical mode has parabolic dispersion 
as k — > owing to rapidly decaying interatomic forces 
for transversal displacements--. Another feature is the 
overlap of the lowest transversal optical mode with the 
acoustical modes. 

In Fig. [3] (b) wc show the phonon dispersion curve of 
BN atomic chain. Two TA modes have low frequency and 
get very small but negative values near the zone center. 
This indicates structural instabilty as A — * oo. However, 
the linear segments of BN atomic chain can be stable. 
Similar to h-BN, the doubly degenerate TO branch over- 
laps with the LA branch. 

For 2D BN honeycomb structure, the unit cell consists 
of two atoms. Accordingly, there are three acoustical 
and three optical branches in Fig. [3] (c). The symme- 
try point group is D3/,. (space group (P-62m)). Optical 
phonon modes at the F-point is given by Aj +2 E . The 
mode A2 is IR active and the E mode is both IR and 
Raman active. The similarity between calculated phonon 
dispersion curves of h-BN and 2D-BN is remarkable. 

We also calculate the phonon dispersion curves of 2D 
BN honeycomb structure by using PAW potentials'^ as 
implemented in VASP'"^ for further checking of the re- 
sults of our phonon calculation. Force constants are de- 
termined from the (8x8x1) supercells. The phonon 
modes were calculated by using the direct method as 
implemented in the PHON^ software. The calculated 
phonon frequencies are almost identical with those cal- 
culated by DFPT method. In Fig. [3] (d), we present 
the phonon density of states calculated for 2D BN hon- 
eycomb structure. Note that both calculations yield that 
TA (or ZA) mode displaying parabolic dispersion gets 
negative frequencies as fe ^ 0. similar to BN atomic 
chains, this indicates structural instability as A — > 00. 
Accordingly, finite size of 2D BN flakes are expected to 
have stable structure. 



V. ID BN NANORIBBONS 

Similar to graphene^, two unique orientation in 2D 
BN yield nanoribbons with uniform edges: These are 
armchair (A-BNNR) and zigzag (Z-BNNR) nanoribbons. 
The profile of the atomic configuration at both edges 
of the nanoribbon determines their electronic and mag- 
netic properties. The properties can be modified by the 
passivation of dangling bond of edge atoms by hydro- 
gen. Because of their interesting electronic and spintronic 
properties, BN nanoribbons are attractive nanostruc- 
tures for various device applications. Electronic proper- 
ties of BN nanoribbons have been investigated in recent 
papersiSiiiii^. Present study is complementary to previ- 
ous studies. 
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FIG. 3: (Color online) Calculated phonon frequencies versus 
k-vectors. (a) h-BN crystal, (b) ID BN atomic chain, (c) 
2D BN honeycomb structure. Phonon modes calculated by 
force constant direct method are shown by the blue-dashed 
curve, (d) Density of phonon frequencies (DOS) for the 2D 
BN honeycomb structure. 



A. Electronic structure 

Here we present the results of our study on the elec- 
tronic and magnetic properties of bare and hydrogen 
passivated A-BNNR and Z-BNNRs. Bare and hydrogen 
passivation A-BNNR are wide band gap semiconductors. 
Similarly, hydrogen passivated Z-BNNRs are also semi- 
conductor. The band gaps of these BN nanoribbons de- 
pend on the width of the nanoribbons w or the numbers 
of BN pairs, n in the primitive unit cell. The variation 
of the band gap Eq as a function of n is given in Fig. U) 
Normally, the properties of nanoribbons approaches to 
those of 2D honeycomb structure as the width n — > cx). 
However, due to the localized edge states the band gap 
of Z-BNNR approaches to a gap smaller than that of 2D 
BN honeycomb structureii. For narrow (n < 8) bare 
and hydrogen passivated A-BNNRs the band gaps vary 
with n, but they are practically unaltered for n > 8. For 
n > 8 the band gap of bare A-BNNR is 0.4 eV smaller 
than that of hydrogen passivated A-BNNR. The band 
gap of hydrogen passivated Z-BNNR is 4.5 eV for n = 3, 
but decrease to 3.8 eV for n = 16. However, its variation 
with n is not monotonic for 5 < n < 13, it rather display 
family dependent oscillatory variation with changes as 
large as 0.4 eV between two consecutive values of n. On 
the other hand, bare Z-BNNRs are found to be metallic. 

The atomic and electronic structure of bare and hy- 
drogen passivated A-BNNR are described in Fig. [5] for 
n — 12. The atoms at the edges of the bare A-BNNR are 
reconstructed; while one edge atom, B is lowering, adja- 
cent edge atom, N is raised. Two bands of edge states 
occur below the conduction band edge. These bands are 
normally degenerate for large n, but split around the 
center of BZ due to their coupling. The bands of edge 
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FIG. 4: (Color online). Energy band gap versus the width of 
the nanoribbons given in terms of the number of B-N atom 
pairs in the primitive unit cell, n. Bare armchair nanorib- 
bons A-BNNR, hydrogen passivated A-BNNR, and hydrogen 
passivated zigzag nanoribbons Z-BNNR. Dotted line indicates 
the bulk band gap of 2D-BN. 




FIG. 5: (Color online) (a) Energy band structure of bare 
armchair nanoribbon A-BNNR having n — 12 B-N pairs in 
the primitive unit cell. At the right hand side of bands, the 
schematic description of atomic structure with primitive unit 
cell delineated by dotted lines and isosurface charge distribu- 
tion of specific states are shown, (b) Same as (a) but the dan- 
gling bonds at both edges are passivated by hydrogen atoms. 



states occur '-^-l eV below the top of the valance band 
edge. Normal states, on the other hand, have charge dis- 
tributed uniformly in the ribbon. Because of the edge 
states the band gap is indirect and is ~4.22 eV wide. 
Upon passivation of the dangling bonds of B and N atoms 
situated at the edges with hydrogen atoms, these edge 
state bands are discarded from the band gap and re- 
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FIG. 6: (Color online) Top panels: Atomic structures of 
zigzag nanoribbons (Z-BNNR). The primitive unit cell has 
n — 6 B-N pairs delineated by dotted lines. The unit cell 
is doubled due to antiferromagnetic interaction between ad- 
jacent N atoms. Middle panels: Isosurface plots of difi'er- 
ence charge density between up spin and down spin states, 
Ap — p(t) — p(i). Bottom panels: Energy band structure 
with dotted (blue) and solid (red) lines showing spin up and 
spin down states, respectively, (a) Bare Z-BNNR; (b) B-side 
free, but N-side is passivated by hydrogen atoms; (c) N-side 
free, but B-side is saturated by hydrogen atoms; (d) Both 
sides are saturated by hydrogen atoms. The bands in (a), 
(b), and (c) are calculated using double cell. 



construction of edge atoms disappear. At the end, the 
band gap of H-passivated A-BNNR becomes direct and 
increases by ~0.3 eV. 

The electronic and magnetic states of Z-BNNR depend 
on whether their edges are passivated with hydrogen 
atoms. While a bare Z-BNNR is magnetic and metallic, 
it becomes nonmagnetic and a wide band gap semicon- 
ductor upon the passivation of B and N atoms at both 
edges. Moreover, its electronic and magnetic properties 
depend on whether only B- or N-side is passivated with 
hydrogen atoms. Accordingly, Z-BNNRs provide us for 
several alternatives for different electronic and magnetic 
properties^*. However, different magnetic states corre- 
sponding to different edge configuration, namely bare or 
hydrogen passivated, are very sensitive to the parame- 
ters of calculation. In Fig. [5] we present the calculated 
electronic structures of a Z-BNNR with n=6 B-N pairs 
in a primitive unit cell for four different cases. These are 
both side free, only N-side is passivated with hydrogen, 
only B-side is passivated with hydrogen and both edges 



are passivated with hydrogen. 

Bare Z-BNNR having both edges are free display differ- 
ent magnetic states (magnetic order), which are close in 
energy. Moreover, the ordering of these magnetic states 
with respect to their energy is sensitive to the criterion 
of energy convergence. To ensure the antiferromagnetic 
(AFM) order at edges, we considered double cells. The 
possible magnetic states are spin-up, spin-down for ad- 
jacent B atoms at one side and spin-up, spin-up for the 
adjacent N atoms at the other side; namely ti / TT spin 
configuration. Other possible spin configurations are ff 

/ ii; TT / TT ; TT / Ti; Ti / TT- We found that the 

spin configuration, fj / TT for Z-BNNR having 12 B- 
N pairs in double unit cell corresponds to the ground 
state. The other excited configurations, TT / ii; TT / TT 
; TT / Ti; Ti / Ti, have 6,7,35,131 meV higher energies 
than ground state. The ordering of these configuration 
is slightly different from that reported earlier^S.. Nev- 
ertheless, the difference between the earlier and present 
ground state energies are within the accuracy limits of 
DFT calculations. The ground state spin configuration 
Ti / TT of the bare Z-BNNR having both edges free is 
found to be ferrimagnetic metal with yu = 1.77 per 
double cell. Whereas the excited magnetic state with 
configuration |T / ii is half metallic. 

In Fig. M (b), Z-BNNR with N-edge passivated with 
hydrogen atoms is an AFM semiconductor. The AFM 
edge state is localized at the B-side. When only the B- 
side is passivated with hydrogen atoms the magnetic edge 
state is, this time, localized at the N-side of the ribbon. 
As seen in Fig. [S] (c) the ground state of Z-BNNR is fer- 
romagnetic with /i —2 /is per double cell. Our calcula- 
tions suggest that the nearest neighbor N-N interaction is 
ferromagnetic, the B-B interaction is antiferromagnetic. 
Finally, the Z-BNNR becomes non magnetic, when the 
atoms at both edges are passivated with hydrogen atoms. 
Earlier, Hwan and Louie^^ studied hydrogen passivated 
A-BNNRs and Z-BNNRs with widths up to 10 nm. Our 
results for hydrogen passivated nanoribbons are in good 
agreement with their results, except that our results for 
zigzag ribbons obtained using GGA as well as LDA ex- 
hibit family dependent oscillations for 5 < n < 13. 



B. Elastic properties 
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FIG. 7; (Color online) (a) Variation of total energy of hydro- 
gen saturated A-BNNR with strain , e is shown by dashed 
curve with large black dots indicating the calculated data 
points (cs = 5c and n — 9). Harmonic, anharmonic and 
plastic regions are distinguished. The harmonic part is fitted 
to a parabola presented by red-solid curve. Atomic structure 
shown by filled, empty and very small empty circles represent 
B, N, and H atoms. Supercell comprising five primitive unit 
cells are shown in the harmonic and plastic regions, (b) Vari- 
ation of K = d^Er/dc^ versus ribbon width n calculated for 
bare and hydrogen passivated A-BNNR. 



The elastic properties of BNNRs are examined through 
the variation of the total energy Et with respect to the 
applied uniaxial strain e = Ac/c^ c being the lattice con- 
stant along the nanoribbon axis. Owing to ambiguities 
in defining the cross section of the ribbon one cannot 
determine the Young's modulus rigorously. Instead we 
calculate k = d'^ET/dc^ from the variation of Et to 
specify the elastic properties of quasi ID nanoribbons. 
In Fig. [7] (a) we show the variation of the total energy 
Et versus e. In order to lift the constraints imposed by 
periodic boundary conditions, calculations are performed 
for a supercell comprising five primitive unit cells having 



lattice constant Cg — 5c. For e < 0.10, the variation 
of Et{() is parabolic, and hence k is independent of e. 
For e > 10 Et{£) curve deviates from parabola and be- 
comes anharmonic. For higher values of strain in the 
plastic region, the ribbon undergoes structural transfor- 
mation. For example, such a transformation occurred 
at e = 0.24 with a sudden change in i?T(e) curve. The 
corresponding structure is illustrated as inset. The lat- 
tice constant Cs increased from the initial value 21.5 A to 
27.4 Acorresponding to e = 0.27. 

In Fig. [7] (b) K versus the width of the ribbon in terms 
of the number of B-N pair in the primitive unit cell n is 
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FIG. 8: (Color online) Variation of total energy of hydro- 
gen saturated Z-BNNR shown by dashed curve with large 
black dots indicating the calculated data points. Harmonic, 
anharmonic and plastic regions are distinguished. The har- 
monic part is fitted to a parabola presented by a red solid 
curve. Atomic structure of the ribbon in a supercell compris- 
ing eight unit cells (cs — 8c and n = 6) are shown before and 
after structural transformation as inset. 



plotted for bare and hydrogen passivated A-BNNR. K{n) 
shows an approximately linear variation indicating that 
the force constant is directly proportional to the width of 
the ribbon. One also sees that the strength of the ribbon 
increases upon passivation with hydrogen. 

The behavior of bare and hydrogen passivated Z- 
BNNR under uniaxial tensile stress is similar to that of 
A-BNNR. In Fig. [5] three regions, namely elastic har- 
monic, elastic- anharmonic and plastic regions are seen. 
The sudden change in the Exie) curve at e 0.23 indi- 
cates a structural phase transformation, where the lat- 
tice constant Cg elongates from the initial e = value of 
19.8 A to 25.7 A corresponding to e = 0.3. The struc- 
ture of hydrogen passivated Z-BNNR before and after 
the structural transformation are shown as inset. Varia- 
tion of K versus the ribbon width n is calculated for bare 
and hydrogen passivated Z-BNNR show an overall linear 
behavior as presented in Fig. [7] (b) . 



C. Vacancy and antisite defects 

It is known that the vacancy defect in 2D 
graphene^Si^i^i^!^ and graphene nanoribbons^i*^ give 
rise to crucial changes in the electronic and magnetic 
structure. According to Lieb's theorem^^, the net mag- 
netic moment per cell is determined with the difference 
in the number of atoms belonging to different sublattices, 
i.e. /i = {Nb — Nis[)fj,B- While DFT calculations on va- 
cancies in 2D graphene and armchair graphene nanorib- 
bons confirmed Lieb's theorem, results are diversified for 
vacancies in zigzag graphene nanoribbons^i*^. There- 
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FIG. 9: (Color online) Relaxed atomic structures and corre- 
sponding energy bands of hydrogen passivated A-BNNR with 
n = 12 having a point defect located periodically in every four 
primitive cell. Blue filled, empty and small circles represent 
B, N, and H atoms, respectively. Blue-dark and yellow-light 
isosurface plots are for spin-up and spin-down states, (a) Sin- 
gle B-vacancy; (b) single N vacancy; (c) B-N divacancy; (d) 
Antisite defect. 



fore, the effect of vacancy defects on the properties of 
BNNRs is of interest. 

Earlier activation energies and reaction paths for diffu- 
sion and nucleation mono and divacancy in h-BN layers 
have been investigated by using density functional tight- 
binding method^^. The formation energies were calcu- 
lated to be 11.22 eV and 8.91 eV, for a B- and N-vacancy, 
respectively. The possible magnetism induced by non- 
magnetic impurities and vacancy defects in a BN sheet 
have been investigated from the first-principles. The 
magnetic moment associated by nonmagnetic atoms sub- 
stituting B or N has been calculated to be l/us^. Based 
on first-principles calculations, the magnetic moment of 
a N-vacancy in a 2D BN sheet has been predicted to be 
1 hb- In the case of a B-vacancy, three neighboring N 
atoms are displaced further apart from each other and 
the net magnetic moment is predicted to be 3 fiB^- An- 
other calculation of defects in a BN monolayer found that 
three dangling bonds associated with a B-vacancy lead to 
total spin S=3/2, i.e 3 /is^''- 

The effects of vacancies of BN nanoribbons have not 
been treated yet. Here we investigated the effect of B-, N- 
, B-|-N-divacancy and B-|-N-anti site on the electronic and 
magnetic properties of A- and Z-BNNR. Within periodic 
boundary conditions, a vacancy defect in an A-BNNR 
of width n = 12 is repeated in every 5th primitive unit 
cell to yield minute defect-defect coupling. As shown in 
Fig. [HI (a) A-BNNR with B-vacancy becomes ferromag- 
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FIG. 10: (Color online) Relaxed atomic structures of hydro- 
gen passivated Z-BNNR with n = 6 having a vacancy defect 
located periodically in every eight primitive cell (cs ~ 8c and 
n = 6). Filled, empty and small circles represent B, N, and H 
atoms, respectively. Blue and pink isosurface plots are spin 
up and spin down states, respectively, (a) Single B-vacancy; 
(b) single N vacancy; (c) B-f N divacancy; (d) anti site defect. 



netic with a net magnetic moment of = 1 hb per unit 
cell. Similarly, a N vacancy gives rise to a net magnetic 
moment oi ^ = 1 per unit cell. A-BNNR having ei- 
ther periodic B-f N-divacancy or anti site defect for every 
five unit cell remain nonmagnetic. The calculated values 
of magnetic moments are in compliance with Lieb's theo- 
rem. We found that the structural relaxation is crucial to 
obtain correct values of magnetic moments. In particu- 
lar, initially we calculated /i = 3 ^b for relaxed structure 
of the B-vacancy. However, the neighboring N atoms 
distorted slightly from their equilibrium, the structure is 
relaxed further and had lowered the total energy. As a 
result, the magnetic moment was calculated as /i = 1 /i^. 
The energy band structures in Fig. [S] (a)-(d) are calcu- 
lated for periodic vacancy defects repeating in every four 
primitive cell. The Fermi levels are assigned according to 
the occupancy of vacancy states. We note that the empty 
state associated with the B-vacancy in Fig. IHlis hole like. 
The states associated with the N- vacancy occur near the 
edge of the conduction band are donor like. 

The situation with Z-BNNR is similar to that in A- 
BNNR discussed above, since hydrogen passivated Z- 
BNNR is nonmagnetic as A-BNNR. A periodic B- or 
N-vacancy repeated in every eight unit cell of hydrogen 
saturated Z-BNNR with n=6 has a net magnetic moment 
of fj, = 1 fiB per supercell. Whereas Z-BNNR passivated 



with hydrogen atoms at both edges and having either 
periodic B-|-N-divacancy or anti site defect repeating in 
every eight unit cell is nonmagnetic. The type of the pe- 
riodic vacancy defect modifies the band gap of Z-BNNR 
from 4.2 eV to 2.21 eV for divacancy, but to 2.8 eV for 
anti site. The calculated magnetic moment of hydrogen 
passivated Z-BNNR are in agreement with Lieb's theo- 
rem. We note that in zigzag graphene nanoribbons mag- 
netic edge states survive even after hydrogen passivation, 
and interact with the magnetic moments of vacancies'^^. 
This interaction causes deviation from the prediction of 
Lieb's theorem. 



VI. DISCUSSION AND CONCLUSIONS 

In various allotropic forms of BN the dimensionality 
play a crucial role. For the sake of comparison, we present 
the calculated values of BN for different allotropic forms 
of BN in different dimensionalities. One sees that the B- 
N double bond of ID BN atomic chain is shortest and is 
1.31 A. The sp2 bond of h-BN and 2D BN has intermedi- 
ate value of 1.45 A. Therefore, h-BN can considered to be 
quasi two dimensional. Three dimensional wurtzite and 
zincblende BN crystal have sp'^-bonding with d — 1.56 A, 
which is largest among the allotropic forms studied here. 
According to GGA results the cohesive energy of 2D BN 
is 3 meV larger than that of h-BN. This is due to fact that 
the GGA calculation cannot account the van der Waals 
interaction between atomic layers of h-BN. However, the 
calculations using LDA, where the van der Waals inter- 
actions are better accounted, yield the cohesive energy 
of h-BN is '^57 meV larger than that of 2D BN as one 
expects. The charge transfer AQ from B to N atom 
increases with decreasing dimensionality. This due to 
fact that d decreases with decreasing dimensionality. As 
for the coordination number increases with increasing di- 
mension. 

In 2D BN honeycomb structures and in its zigzag and 
armchair nanoribbons, the B-N bond formed from the 
bonding sp^ hybrid orbitals from B and N atoms is es- 
sential. Owing to the transfer of charge from B to N the 
B-N bond acquires an ionic character, which underlies 
the semiconducting properties with wide band gap. 

Bare armchair nanoribbon of 2D BN is again a non- 
magnetic wide band gap semiconductor, the band gap of 
which is practically unaltered with width n > 8. Upon 
passivation with hydrogen band gap of the ribbon in- 
crease by 0.3 eV. As for zigzag nanoribbons, they provide 
a number of interesting properties. When its both edges 
are bare, it is ferromagnetic metal. When its N-edge is 
passivated with hydrogen, it becomes an antiferromag- 
netic semiconductor. In the reverse case, namely when 
B-side is passivated, it becomes a ferromagnetic semicon- 
ductor. When both edges are passivated, it becomes a 
nonmagnetic, wide band gap semiconductor. The band 
gap as well as the magnetic state of a ribbon can be mod- 
ified by periodic vacancy defects. Finally, BN nanorib- 



10 



TABLE I: Values of the bond length d in A, cohesive energy 
Ec in eV per B-N pair, band gap Eg in eV, charge transfer 
from B to N (AQ) in electrons and lattice constants (a,c) 
in Acalculated for various allotropic forms of BN in different 
dimensionality. 





d 


Ec 


Eg 


A Q 


Lattice 


ID Chain 


1.307 


-16.04 


3.99 


0.511 


a=2.614 


2D BN 


1.452 


-17.65 


4.64 


0.429 


a=2.511 


h-BN 


1.450 


-17.65 


4.47 


0.416 


a=2.511 , c/a=2.66 


Wurtzite 


1.561 


-17.45 


5.726 


0.342 


a=2.542 , c/a=1.63 


Zincblende 


1.568 


-17.49 


4.50 


0.334 


a=2.561 



bons have been found to be strong, quasi one dimen- 
sional and stable structures. They can sustain up to high 
strains, and they stretch in the plastic region with struc- 
tural transformations. 

Briefly, the calculated electronic, magnetic and me- 
chanical properties of 2D BN honeycomb structure and 



its nanoribbons present interesting but some differences 
from graphene. In this respect BN honeycomb structure 
and its nanoribbons are complimentary to graphene. The 
properties of 2D BN honeycomb structure can be changed 
upon functionalization with foreign atoms. Interesting 
quantum structures (such as single and series quantum 
dots, resonant tunneling double barriers and multiple 
quantum well structures) based on heterostructures and 
core shell structures of lattice matched graphene and BN 
can be formed, since the band gap of BN is much larger 
than that of graphene. 
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